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I.  INTRODUCTION 


The  ability  to  compute  the  base  region  flow  field  for  projectile  configurations  using 
Navier-Stokes  computational  techniques  has  been  developed  over  the  past  few  yearsltJi3 
This  capability  is  very  important  for  determining  aerodynamic  coefficient  data  including 
the  total  aerodynamic  drag.  The  majority  of  base  flow  calculations  to  date  have  modeled 
the  base  region  as  a  flat  solid  surface.  Many  of  the  actual  configurations  have  some  form  of 
base  cavity.  General  opinion  has  been  that  the  inclusion  of  a  base  cavity  or  modifications 
to  the  interior  cavity  of  a  projectile  base  would  have  little  or  no  effect  on  the  overall  flight 
performance  parameters. 

The  M825  projectile  under  certain  conditions  is  expected  to  be  acroballistically  similar 
to  its  parent  configuration  the  M4S3A1.  The  M825  has  an  aluminum/steel  base  which  is 
configured  as  a  flat  cavity  (standard).  A  recent  Product  Improvement  Program  (PIP), 
undertaken  to  reduce  the  production  costs  and  improve  shell  integrity,  resulted  in  the 
design  of  a  new  base  configuration.  This  new  PIP  configuration  has  an  all  Bteel  base  and 
contains  a  dome  cavity.  A  series  of  aeroballistic  tests  4  were  conducted  in  the  Transonic 
Range  Facility  of  the  US  Army  Ballistic  Research  Laboratory  (BRL),  to  determine  any 
difference  in  the  aeroballistics  which  may  occur  between  the  standard  and  dome  base 
configurations.  As  a  result  of  these  tests,  differences  in  aerodynamic  performance  were 
found  to  exist  between  the  two  rounds.  The  most  significant  changes  in  the  aerodynamic 
data  were  in  the  lift  and  static  moment  coefficients.  The  drag  was  found  to  differ  by  a  few 
percent  with  the  dome  configuration  having  the  lower  drag  at  low  transonic  speeds. 

A  computational  study  was  undertaken  to  determine  the  ability  of  the  present  Navier- 
Stokes  codes  to  predict  these  differences  and  to  further  understand  the  fluid  dynamic 
behavior  which  can  account  for  such  small  changes.  The  use  of  Navier-Stokes  codes  can 
provide  a  detailed  description  of  the  flow  field  associated  with  the  M825  configuration 
as  well  as  the  integrated  aerodynamic  coefficients.  The  initial  results  for  the  zero  degree 
angle  of  attack  case  have  been  reported  by  Sahu  et.  al.5  .  This  was  accomplished  using  an 
axisymmetric  base  flow  code  and  the  results  showed  the  same  effect  as  the  range  data,  that 
is,  a  small  reduction  in  the  total  aerodynamic  drag  at  low  transonic  speeds  (  M  <  0.95) 
for  the  dome  base  configuration.  The  trend  reversed  at  high  transonic  speeds  (  M  >  0.98). 
This  report  describes  an  extension  of  that  work  into  three  dimensions. 

Numerical  computations  have  been  performed  using  a  3D  zonal,  implicit  Navier-Stokes 
code.  The  Mach  number  range  was  0.8  <  M  <  1.5  for  an  angle  of  attack,  a  =  4.0°  . 
Results  presented  include  the  qualitative  features  of  the  base  region  flow  field  for  the  two 
base  cavities.  Aerodynamic  force  and  moment  coefficients  have  been  obtained  from  the 
computed  solutions  and  are  presented  as  a  function  of  Mach  number.  Computed  results 
show  small  differences  in  normal  force  and  pitching  moment  coefficients  similar  to  that 
found  in  the  range  data. 
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II.  GOVERN -NG  EQUATIONS  AND  SOLUTION 

TECHNIQUE 


The  complete  set  of  time-dependent  thin-layer  Navier-Stokes  equations  is  solved  nu¬ 
merically  to  obtain  a  solution  to  this  problem.  The  numerical  technique  used  is  an  implicit 
finite  difference  scheme.  Although  time-dependent  calculations  are  made,  the  transient 
flow  is  not  of  primary  interest  at  the  present  time.  The  steady  flow,  which  is  the  desired 
result,  is  obtained  in  a  time  asymptotic  fashion. 


1.  GOVERNING  EQUATIONS 


The  complete  set  of  three  dimensional,  time  dependent,  generalized  geometry,  thin- 
layer,  Navier-Stokes  equations  for  general  spatial  coordinates  £,  r?,  (  can  be  written  as6: 


where 


and 


drq  +  d,F  +  <9„G  +  dcH  =  Rt~'dcS 
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and  where 

0 

MCx  +  C  +  C*W  +  ^(G«<  +  Cyyc  +  G“><)G 
MCx  +  Cy  +  C  z)V(  +  ^(G«<  +  Cy^C  +  C*^<)Cy 
M(Cx  +  Cv  +  Cx)WC  +  ^(G^C  +  CyV(  +  Gt°<)G 
{(Cx  +  CI  +  Cx)[|(u2  +  v2  +  W7)c 

+ - ifi - ] 

Pr(7-l)J 

+  ^(G«  +  Cyt>  +  C*w)(C*Uc  +  CyU(  +  G^c)} 

In  equation  (1),  the  thin-layer  approximation  is  used  and  the  viscous  terms  involving 
velocity  gradients  in  both  the  longitudinal  and  circumferential  directions  are  neglected. 
The  viscous  terms  are  retained,  however,  for  velocity  gradients  in  a  direction  nearly  normal 
to  the  surface  where  large  fiowfield  gradients  exist.  These  viscous  terms  in  £  are  collected 
into  the  vector  S. 

For  this  computation,  the  diffusion  coefficients  p  and  k  contain  molecular  and  turbu¬ 
lent  parts.  The  turbulent  contributions  are  supplied  through  an  algebraic  eddy-viscosity 
hypothesis  which  has  been  developed  by  Baldwin  and  Lomax.7 

The  velocities  in  the  £,  77,  and  £  coordinate  directions  can  be  written 

U  =  £*  +  «£*  +  +  u>G 

V  =  r/t  +  urjx  +  vt)v  +  wrix 

W  =  G  +  "G  +  t>Cy  +  wCz 

which  represent  the  contravariant  velocity  components. 

The  Cartesian  velocity  components  (u,  v,  w)  are  retained  as  the  dependent  variables 

and  are  nondimensionalized  with  respect  to  (the  free  stream  speed  of  sound).  The  local 

pressure  is  determined  using  the  relation 

P  =  (7  —  l)[c  -  0.5p(u2  +  v7  4-  u>2)]  (4) 

where  7  is  the  ratio  of  specific  heats.  Density  ( p )  is  referenced  to  Poo  and  the  total  energy 
(e)  to  poo^l o*  The  transport  coefficients  are  also  nondimensionalized  with  respect  to  the 
corresponding  free  stream  variables.  Thus  the  Prandtl  number  which  appears  in  S  is 
defined  as  Pr  =  cp00p0 o/k^. 

In  differencing  these  equations  it  is  often  advantageous  to  difference  about  a  known 
base  solution  denoted  by  subscript  0  as 
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(5) 


$t{Q  —  Q  o)  +<5{(F  —  Fo)  +  6V(G  —  Go)  +  &c(H  —  Ho)  —  Re~16((S  —  So) 

—  ~dTQo  —  d^Fo  ~  5,, Go  —  d(Ho  +  Re~id(So 

where  6  indicates  a  general  difference  operator,  and  d  is  the  differential  operator.  If  the 
base  state  can  be  properly  chosen,  the  differenced  quantities  can  have  smaller  and  smoother 
variation  and  therefore  less  differencing  error. 


2.  NUMERICAL  TECHNIQUE 

The  implicit  approximately  factored  scheme  for  the  thin  layer  Navier  Stokes  equations 
that  uses  central  differencing  in  the  tj  and  (  directions  and  upwinding  in  £  is  written  in 
the  form 


[i  +  M{(A+)r  +  h6(Cn  -  hRt~%J~lMnJ  -  Ale] 

x  [/  +  h6{ (A- )n  +  hS.B”  -  A|n]  A Qn  =  (6) 

-At{*;[(F+r  -  F+)  +  6[[(F-)n  -  F-]  +  *„(Gn  -  Goo) 

+*<(#"  -  Hoo)  -  Fe",^(Sn  -  Soo)}  -  De(Qn  -  QTO) 


where  h  =  At  or  (At  )/2  and  the  free  stream  base  solution  is  used.  Here  S  is  typically  a 
three  point  second  order  accurate  central  difference  operator,  6  is  a  midpoint  operator  used 
with  the  viscous  terms,  and  the  operators  and  6^  are  backward  and  forward  three-point 

A  A  *  A  A 

difference  operators.  The  flux  F  has  been  eigensplit  and  the  matrices  A,  B,C,  and  M 
result  from  local  linearization  of  the  fluxes  about  the  previous  time  level.  Here  J  denotes 
the  Jacobian  of  the  coordinate  transformation.  Dissipation  operators,  Dt  and  A  are  used 
in  the  central  space  differencing  directions. 

The  smoothing  terms  used  in  the  present  study  are  of  the  form: 

•o.i,  =  p{B)0s + 

.A  I,  =  (At)J~l[e26p(B)P6  +  2.5^6 p{B)S]J 

where  0  =  and  where  p{B)  is  the  true  spectral  radius  of  B.  The  idea  here  is  that 

the  fourth  difference  will  be  tuned  down  near  shocks,  that  is,  as  0  gets  large  the  weight  on 
the  fourth  difference  drops  down  while  the  second  difference  times  up. 

For  simplicity,  all  the  boundary  conditions  have  been  imposed  explicitly.  On  the  body 
surface,  the  no-slip  boundary  condition  is  used  and  the  wall  temperature  is  specified.  Free 
stream  boundary  conditions  are  used  at  the  computational  outer  boundary.  A  symmetry 


4 


boundary  condition  is  imposed  at  the  circumferential  edges  of  the  grid  while  a  simple 
extrapolation  is  used  at  the  downstream  boundary.  The  flowfield  is  initially  set  to  free 
stream  conditions  everywhere  and  then  advanced  in  time  until  a  steady  state  solution  is 
obtained.  Atmospheric  flight  conditions  were  used. 


3.  COMPOSITE  GRID  SCHEME 

In  the  present  work,  a  simple  composite  grid  scheme8  has  been  used  where  a  large 
single  grid  is  split  into  a  number  of  smaller  grids  so  that  computations  can  be  performed 
on  each  of  these  grids  separately.  These  grids  use  the  available  core  memory  one  grid 
at  a  time,  while  the  remaining  grids  are  stored  on  an  external  disk  storage  device  such 
as  the  solid  state  disk  device  (SSD)  of  the  Cray  X-MP/48  computer.  The  Cray-2  has  a 
large  incore  memory  to  fit  the  large  single  grid.  However,  for  accurate  geometric  modeling 
of  complex  projectile  configurations  which  include  blunt  noses,  sharp  comers  and  base 
cavities,  it  is  also  desirable  to  split  the  large  data  base  into  a  few  smaller  zones  on  Cray-2 
as  well. 

The  use  of  a  composite  grid  scheme  requires  special  care  in  storing  and  fetching  the 
interface  boundary  data  ,  i.e.,  the  communication  between  the  various  zones.  In  the  present 
scheme,  there  is  a  one  to  one  mapping  of  the  grid  points  at  the  interface  boundaries  and 
thus,  no  interpolations  are  required.  Details  of  the  data  storage,  data  transfer  and  other 
pertinent  information  such  as  metric  and  differencing  accuracy  at  the  interfaces  can  be 
found  in  Reference  S  and  9.  This  scheme  has  been  successfully  used  by  Sahu9  to  compute 
three  dimensional  transonic  flow  over  two  projectiles  .  The  computed  results  clearly  showed 
the  transonic  critical  aerodynamic  behavior  in  pitching  moment  coefficient  observed  in  free 
flights.  The  present  work  is  a  further  application  of  this  technique  to  a  more  complicated 
projectile  with  base  cavities. 
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III.  MODEL  GEOMETRY  AND  COMPUTATIONAL  GRID 


The  external  configuration  of  the  M825,  excluding  the  base,  is  similar  to  the  M483A1 
shown  in  Figure  1.  The  features  of  this  projectile  which  have  not  been  modeled  exactly 
are  the  meplat  on  the  fuze  and  the  rotating  band  near  the  base.  The  rotating  band 
was  eliminated  for  simplicity  and  the  meplat  was  modeled  as  a  hemisphere  cap.  The 
computational  model  is  shown  in  Figure  2  and  consists  of  a  2.84  caliber  nose,  a  2.7  caliber 
cylindrical  section,  and  a  0.26  caliber  8°  boattail.  The  ogive  contour  as  well  as  the  undercut 
on  the  cylindrical  section  were  matched. 

The  current  problem  of  interest  is  the  effect  of  the  different  base  geometries  on  the 
overall  projectile  aerodynamics.  Figure  3  shows  the  standard  and  dome  base  configurations. 
The  standard  base  is  a  combination  of  aluminum  and  steel  and  contains  a  base  cavity 
which  is  characterized  as  a  flat  surface.  The  PIP  configuration  is  an  all  steel  base  and 
is  characterized  as  a  dome  surface.  The  cavity  volume  is  also  significantly  larger  for  the 
dome  configuration. 

The  solution  technique  requires  the  discretization  of  the  entire  flow  region  of  interest 
into  a  suitable  computational  grid.  The  grid  outer  boundary  has  been  placed  at  2.5  body 
lengths  upstream  and  surrounding  the  projectile.  The  downstream  boundary  was  placed 
at  2  body  lengths.  Since  the  calculations  are  in  the  subsonic/transonic  regime  the  compu¬ 
tational  boundaries  must  extend  out  beyond  the  influence  of  the  body.  This  ensures  that 
the  boundary  conditions  specified  in  the  flow  code  are  satisfied. 

Figure  4  and  Figure  5  show  the  grids  generated  for  the  standard  base  and  dome  base 
configurations,  respectively.  Each  of  these  grids  consists  of  225  points  in  the  streamwise 
direction  and  50  points  in  the  normal  direction.  This  is  broken  down  into  two  sections: 
a  body  region  and  a  base  region.  The  surface  points  for  each  region  are  selected  using 
an  interactive  design  program.  Each  grid  section  is  then  computed  separately  using  a 
hyperbolic  grid  generation  program  10 .  Longitudinally,  there  are  106  points  along  the 
projectile  surface  and  60  points  in  the  base  region  downstream  of  the  base  corner.  The 
normal  distribution  of  points  in  base  region  consists  of  50  points  along  the  base  cavity.  An 
expanded  three  dimensional  view  of  the  base  grid  is  shown  in  Figure  .6.  This  grid  has  33 
points  in  the  circumferential  direction.  The  generally  flat  sections  on  the  standard  base 
enabled  a  grid  to  be  routinely  generated.  However,  due  to  the  extreme  concavity  the  grid 
for  the  dome  base  (Figure  7)  required  an  increase  in  the  smoothing  values  used  by  the 
hyperbolic  grid  generator,  as  well  as  the  addition  of  a  grid  averaging  technique. 


IV.  RESULTS 

Numerical  computations  have  been  made  for  both  the  standard  and  the  dome  base 
configurations  for  a  range  of  Mach  numbers  from  M  =  0.80  to  1.5  and  at  4  degrees  angle 
of  attack.  Computed  results  obtained  at  zero  degree  angle  of  attack  are  also  included  for 
comparison  purposes. 

A  few  qualitative  results  are  presented  next.  Figure  8  and  Figure  9  show  the  velocity 
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vectors  in  the  base  region  for  both  base  configurations  at  M^O.98  and  a  =  0.0°.  The 
recirculatory  flow  in  the  base  region  is  evident  and  as  expected,  is  symmetric.  As  shown 
in  Figure  8,  the  recirculation  region  for  the  standard  base  extends  to  about  one  and  a  half 
caliber  downstream  of  the  base  corner.  The  back  flow,  upon  reaching  the  cavity  follows 
the  contour  of  the  cavity  and  leaves  the  cavity  pushing  the  flow  upwards.  The  shear  layer 
leaving  the  base  corner  is  displaced  upwards  weakening  the  expansion  at  the  base.  Figure  9 
for  the  dome  configuration  shows  a  weak  secondary  bubble  inside  the  cavity  in  addition 
to  the  primary  bubble.  The  flow  again  follows  the  contour  of  the  cavity  and,  upon  leaving 
the  dome  cavity,  is  almost  parallel  to  the  streamwise  direction.  This  flow,  thus,  has  less 
effect  on  the  free  shear  layer  and  doesn’t  weaken  the  expansion  at  the  base  corner  as  much 
compared  to  the  standard  base.  The  net  effect  is  that  the  size  of  the  primary  bubble  for 
the  dome  base  is  slightly  smaller  than  that  for  the  standard  base.  The  reattachment  point 
is  therefore  closer  to  the  base  and  results  in  lower  base  pressure  or  higher  base  drag  at 
this  Mach  number5  .  Figure  10  and  Figure  11  show  the  velocity  vectors  in  the  base  region 
for  the  base  configurations  at  11/00=0.98  and  a  =  4.0°  for  both  windside  (bottom  half) 
and  leeside  (top  half).  Again  the  recirculatory  flow  in  the  base  region  is  evident  and  as 
expected,  the  flow  in  the  wake  is  asymmetric.  As  shown  in  these  figures,  the  separation 
bubbles  on  windside  and  leeside  in  the  wake  differ  in  size  and  shape  (  the  one  on  windside 
being  more  thin  and  elongated  ).  In  addition,  a  number  of  secondary  separation  bubbles 
can  be  seen  to  form  inside  the  cavity  for  both  base  configurations.  These  changes  in  the  flow 
structure  contribute  to  small  changes  in  the  base  pressure  and  thus,  to  the  aerodynamic 
forces  and  moments. 

Figure  12  and  Figure  13  show  the  Mach  number  contours  in  the  base  region  for  both 
base  configurations  at  Moc,= 0.98  and  a  =  0.0°.  These  figures  show  the  flow  expansion 
at  the  ogive  corner,  boattail  corner  and  the  base  comer.  One  can  also  see  a  shock  wave 
on  the  cylinder  portion  of  the  projectile  as  well  as  a  recompression  shock  system  which 
exists  downstream  of  the  base  corner.  The  flow  field  is  symmetric  for  this  condition.  As 
angle  of  attack  is  increased  to  4  degrees,  the  flow  field  becomes  asymmetric  (  see  Figure  14 
and  Figure  15  ).  A  small  asymmetry  can  be  observed  in  the  location  of  the  shock  wave 
on  the  cylinder.  The  windside  shock  is  further  aft  compared  to  the  corresponding  one  on 
the  leeside.  The  asymmetry  can  be  clearly  seen  in  the  wake  flow  and  its  associated  shock 
system.  As  can  be  seen  in  these  figures,  the  wake  flow  field  changes  for  the  different  base 
configurations. 

The  entire  flowfield  over  the  projectile  including  the  base  region  is  computed.  There¬ 
fore,  the  computed  results  include  any  upstream  influence  the  base  region  flow  may  have 
on  the  boattail  flowfield.  Surface  pressures  including  the  base  pressure  and  the  viscous 
stresses  are  known  from  the  computed  flow  field  and  can  be  integrated  to  give  the  aerody¬ 
namic  forces  and  moments.  Figure  16  shows  the  computed  base  pressure  distribution  for 
the  dome  base  configuration  at  Af0O=l.l  and  a  —  4.0°  for  windside  and  leeside.  As  seen 
in  this  figure,  the  pressure  on  the  windside  (  Z/D=-Q.b  )  is  higher  than  the  pressure  on 
leeside.  Since  this  pressure  acts  normal  to  the  inside  surface  of  the  cavity,  it  produces  a 
downward  force  (see  Figure  17  ).  Figure  17  shows  the  normal  force  coefficient  for  the  dome 
base  as  a  function  of  Mach  number.  The  dotted  line  represents  the  normal  force  coefficient, 
Cn  for  the  dome  base  projectile  where  the  base  region  is  excluded  in  the  force  and  moment 
calculations.  The  solid  line  is  for  the  entire  projectile  including  the  contribution  from  the 
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base  region.  As  mentioned  earlier,  the  base  region  produces  a  negative  contribution  and 
thus,  the  normal  force  is  reduced  somewhat  for  high  transonic  speeds  (  M  >  0.90  ).  The 
reverse  is  true  for  low  transonic  speeds  (  M  <  0.85  )  .  Figure  18  shows  the  normal  force 
coefficient  comparison  for  both  base  configurations  as  a  function  of  Mach  number.  The 
dotted  line  represents  the  dome  base  result  whereas  the  solid  line  shows  the  result  for  the 
standard  base.  As  seen  in  this  figure,  the  dome  base  has  a  higher  normal  force  coefficient 
at  transonic  speeds  compared  to  the  standard  base.  The  difference  is  small,  of  the  order 
of  a  few  percent,  at  higher  transonic  speeds  (  M  >  0.90  )  and  gets  as  large  as  10-12  %  at 
low  transonic  speeds  (  M  <  0.90  ). 

An  aerodynamic  coefficient  which  is  of  primary  concern  is  the  pitching  moment  co¬ 
efficient,  Cmo.  Figure  19  shows  the  Cma  comparison  for  both  base  configurations.  The 
computed  Cmo  is  also  compared  with  the  range  data  4  for  both  base  configurations.  Here 
Cma  is  referenced  to  the  center  of  gravity  of  the  projectile.  The  computed  result  clearly 
shows  a  sharp  rise  in  Cma  between  M  =  0.80  to  0.88  which  is  followed  by  a  sharp  drop  as 
Mach  number  is  increased  to  M  =  0.95  .  As  the  Mach  number  is  increased  further  Cmo 
rises  gradually  again  which  is  unlike  the  behavior  of  other  projectiles  such  as  the  M549 
projectile  9.  This  critical  behavior  in  Cma  observed  in  the  data  is  clearly  predicted  in  the 
numerical  computations.  The  overall  comparison  of  the  computed  result  with  the  range 
data  is  fair.  As  seen  in  the  range  data,  the  dome  base  configuration  has  lower  Cma  at 
transonic  speeds  compared  to  the  standard  base  configuration  and  this  trend  is  also  seen 
in  the  computed  results. 


V.  CONCLUSIONS 


Three  dimensional  numerical  computations  have  been  made  for  a  projectile  with  two 
base  cavity  configurations  at  transonic  speeds.  Computed  results  show  differences  in  the 
qualitative  features  of  the  base  region  flowfield  between  the  two  base  cavities.  Changes 
in  the  base  cavity  configuration  have  been  found  to  affect  the  normal  force  and  pitching 
moment  coefficient.  Differences  in  these  coefficients  of  between  0  to  12%  have  been  pre¬ 
dicted  and  are  compared  with  the  range  data.  The  dome  base  configuration  produces  a 
higher  normal  force  and  a  lower  pitching  moment  them  the  standard  base  at  these  transonic 
speeds  and  shows  the  same  trend  observed  in  the  range  data. 
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Figure  1:  Actual  XM825  projectile. 


figure  -1  Computational  grid  for  the  standard  base. 


Figure  5:  Computational  grid  for  the  dome  base. 
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Figure  G  Base  region  grid  for  tlie  standard  base 


Figure  7:  Base  region  grid  for  the  dome  base. 
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iguro  12.  Mach  number  contours, .\/K  =  0-9$,  a  =  O.O5,  (standard  base) 
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Figure  1-F  Mad)  number  contours. .\fK  =0.9$,  a  —  \  O'  ,  (standard  base) 


Figure  15:  Mach  number  contours, M&  =0.98,  a  =  4  0°,  (dome  base). 
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Figure  16.  Base  pressure  distribution  for  the  dome  base,Moe,=l.l,  a  =  4.0° 
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Figure  18.  Normal  force  coefficient, CN  vs  Mach  number,  a  =  4.0°,  (standard  base  and 
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6: 'General  Comments  wfcat*do  you  think  should  be  changed  to  improve  future  reports?  (Indicate 
changes,  to  organization,  technical  content,  format,  etc.)  _ _ 


Name 

CURRENT  Organization 

ADDRESS  _ 

Address 

City,  State,  Zip  Code 

7.  If  ind  .eating  a  Change  of  Address  or  Address  Correction,  please  provide  the  New  or  Correct 
Address  in  Block  6  above  and  the  Old  or  Incorrect  address  below. 


Name 

OLD  Organization 

ADDRESS  _ 

Address 

City,  State,  Zip  Code  , 


(Remove  this  sheet,  fold  as  indicated,  staple  or  tape  closed,  and  mail.) 


